Dynamic aspect of the chiral phase transition 
in the mode couphng theory 



K. Ohnishi'^'^, K. Fukushima ^'^ and K. Ohta^ 

'^Institute of Physics, University of Tokyo, 
3-8-1 Komaba, Meguro-ku, Tokyo 153-8902, Japan 

^ Center for Theoretical Physics, Massachusetts Institute of Technology, 
77 Massachusetts Avenue, Cambridge, MA 02139, U.S.A. 

'^Department of Physics, University of Tokyo, 
7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan 



Abstract 

We analyze the dynamic aspect of the chiral phase transition. We apply the mode 
coupling theory to the linear sigma model and derive the kinetic equation for the 
chiral phase transition. We challenge Hohenberg and Halperin's classification scheme 
of dynamic critical phenomena in which the dynamic universality class of the chiral 
phase transition has been identified with that of the antiferromagnet. We point out 
a crucial difference between the chiral dynamics and the antiferromagnet system. 
We also calculate the dynamic critical exponent for the chiral phase transition. 
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1 Introduction 



Among the central issues in the study of QCD (Quantum ChromoDynamics) , 
there exists the problem of the spontaneous chiral symmetry breaking and 
its restoration, i.e., the chiral phase transition. It is believed that the chiral 
phase transition occurs under an extreme condition like high temperature or 
high density. The critical dynamics of the second order chiral phase transi- 
tion for the two flavor massless QCD is not only of theoretical interest but 
also of importance in understanding the relativistic heavy ion collision ex- 
periments. The critical phenomena significantly affect how hot matter cools 
down. There are a lot of works aimed to look into the static (equilibrium) 
nature of the chiral phase transition, including the lattice QCD simulation 
[1,2,3,4,5]. Our objective in the present paper is to investigate the dynamic 
(non-equilibrium or real-time) property of the second order chiral phase tran- 
sition [6,7,8,9,10,11,12,13]. Knowledge of the dynamic critical behavior would 
be indispensable to understand the relativistic heavy ion collision experiments 
because the system evolves essentially through non-equilibrium states. Anoma- 
lous dynamic critical phenomena such as the critical slowing down and the 
softening of propagating modes have significant relevance to the matter going 
toward or coming from a quark-gluon plasma. Throughout this paper, we will 
restrict ourselves to the chiral phase transition at zero baryon number density. 

In general, whenever we consider a critical point of second order, it is useful 
to classify the critical point according to the universality class. The critical 
points in the same universality class are to show exactly the same critical be- 
havior unless there appears another fixed point or the interaction is infinitely 
long-ranged. For the static case, as is well known, the universality class is de- 
termined solely by the global symmetry of the system. The systems having the 
same symmetry belong to the same static universality class. On the other hand, 
the dynamic universality class is more complicated to classify. A classification 
method, which we will briefiy review in the next section, has been proposed 
and established by Hohcnbcrg and Halperin. Actually, all the observed critical 
points, such as the liquid-gas transition, transitions in a ferromagnet, in an 
antiferromagnet and in a superfluid, have been classified systematically and 
comprehensively in the review article by Hohenberg and Halperin [14] . 

Consideration on the universality class of the chiral phase transition has al- 
ready been given. The static universality class of the two-flavor chiral phase 
transition is identifled with that of the ferromagnet and antiferromagnet be- 
cause they all have the rotational 0(4) symmetry [1]. The dynamic universal- 
ity class of the chiral phase transition has also been discussed by Rajagopal 
and Wilczek [7]. They argue that the chiral phase transition belongs to the 
same dynamic universality class as that of the antiferromagnet. Certainly it is 
concluded that the dynamic universality classes of the two systems would be 
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identical according to the classification method of Hohenbcrg and Halpcrin. 
Since the nature of the antiferromagnet has been analyzed and understood 
enough, the identification of the dynamic universality class means that the 
nature of the chiral phase transition is understood as well. There does exist, 
however, a crucial difference between the two systems in the appearance of 
physical modes, as pointed out in Ref. [10]. The difference is so crucial that 
the dynamic universality classes of the two systems cannot be regarded as 
identical any longer and Hohenberg and Halperin's classification method itself 
must be reexamined carefully. In the next section, we will take a close look at 
the two systems to reveal the difference. 

Once the two systems are realized to be distinctive, it is necessary to reanalyze 
the dynamic aspect of the chiral phase transition without resorting to Hohen- 
berg and Halperin's classification method. We will adopt the mode coupling 
theory [15,16,17,18] in order to disclose the physical modes and write down 
the kinetic equations in a chiral effective model. 

The plan of this paper is as follows. In Sec. 2, we will give a general review 
on the dynamic critical phenomena and point out the difference between the 
antiferromagnet and the chiral phase transition. In Sec. 3, we will review the 
mode coupling theory briefly. In Sec. 4, we will apply the mode coupling theory 
to the 0(2) linear sigma model in order to see what kinds of slow modes appear 
near the phase transition. In Sec. 5 we will derive the kinetic equation in the 
0(4) linear sigma model and we will calculate the dynamic critical exponents 
in Sec. 6. The exponent will be found to be 2; = 1 — ry/2 = 0.98, which is 
contrasted with z ~ d/2 of the antiferromagnet. Finally, we will give the 
summary in Sec. 7. 



2 Dynamic universality class of the chiral phase transition 



In this section, we will discuss the dynamic universality class of the chiral 
phase transition with emphasis on the difference from the antiferromagnet. 
For that purpose, it is inevitably necessary to explain the notion of the slow 
mode. In Subsec. 2.1, we will give a general review on the dynamic critical 
phenomena, in which we will explain the slow variable, the slow mode and 
the classification method of the dynamic universality class by Hohenberg and 
Halperin. In Subsec. 2.2, we will discuss the chiral phase transition. 
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2. 1 Dynamic critical phenomena 



One of the prominent ingredients of the dynamic critical phenomena would 
be the critical slowing down, that is, the long (and infinite) relaxation time 
near (and at) a critical point. The critical slowing down is attributed to the 
slow motion of the long wavelength fluctuations of the slow variables. The 
slow variables near a critical point consist of the order parameter and the con- 
served quantities involved in the system. The slowness of the order parameter 
fluctuations can be understood by considering a flat thermodynamic potential 
with respect to the order parameter near the critical point, which indicates 
a weak restoring force on the order parameter fluctuations. Alternatively we 
can imagine the correlated domains of the order parameter to understand the 
critical slowing down. Near the critical point, the size of the domains, i.e., the 
correlation length ^, becomes so large that the motion of the domains gets 
very dull and slow. Since ^ diverges just at the critical point, the slowness of 
the motion becomes infinite, leading to the divergent relaxation time. 

The motion of the long wavelength fluctuation of a conserved quantity is 
kinematically slow even away from the critical point. In fact, the fluctuation 
of a conserved quantity with zero wavenumber, ^=0, has an infinite relaxation 
time because of conservation of the quantity. [Note that d^j_j^ = means j'^ 
is time independent when q = 0.] The fluctuation with a small but flnite 
wavenumber has a long but flnite relaxation time. Near the critical point, it 
is known that the coupling between the order parameter and the conserved 
quantities becomes very strong. Through the coupling, the fluctuations of the 
conserved quantity even with a flnite wavenumber come to have an inflnite 
relaxation time at the critical point. 

We have seen the slowness of the long wavelength fluctuation of the slow vari- 
ables. The relaxation time for those fluctuations is so long and amounts to the 
macroscopic scale belonging to the kinetic regime. On the other hand, for the 
degrees of freedom other than the long wavelength fluctuation of the slow vari- 
ables, namely the short wavelength fluctuations of the slow variables and the 
fluctuations of the coordinates other than the slow variables, the motions re- 
main rapid near the critical point. Consider a non-equilibrium system near the 
critical point after a macroscopic time has elapsed. The non-equilibrium fluc- 
tuations of the rapid degrees of freedom will have already relaxed to constitute 
the heat bath, while the long wavelength fluctuations of the slow variables will 
still remain as non-equilibrium fluctuations. Thus, by tracing the time evolu- 
tion of the long wavelength fluctuations of the slow variables, we can make the 
sufficient and even complete description of non-equilibrium processes in the ki- 
netic regime. The mode coupling theory developed by Kawasaki [15,16,17,18] 
is one of the schemes which provides us with the equation of motion for the 
critical fluctuations, that is, the kinetic equation near the critical point. 
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The long wavelength flTictuations of the slow variables or their linear com- 
binations constitute what is called the slow modes. The slow modes appear 
as narrow peaks in the spectral functions N{uj, q) associated with the slow 
variables. There are two kinds of slow modes; a diffusive mode and a prop- 
agating mode. The diffusive mode is a mode that is purely dissipative, while 
the propagating mode involves an oscillation besides a dissipation. Examples 
for the propagating mode are the soimd wave, the spin wave and the "par- 
ticle mode." For definiteness, consider the critical fluctuation with a fixed 
wavenumber q [q not necessarily vanishing] in the following explanation. 
The diffusive mode corresponds to the pole of the spectral function in the 
complex plane of the frequency uj having only an imaginary part F. In the 
spectrum as a function of c<j, the diffusive mode is seen as a narrow peak with 
the width F, the peak position of which is located at a; = (so called the 
central peak). The width F represents the decay rate of the fluctuation, that 
is, the inverse relaxation time. As the critical point is approached, F goes to 
zero, implementing the critical slowing down. 

As for the propagating mode, the spectral function has the form that the pole 
with respect to complex uj has the real part Q as well as the imaginary part F. 
The peak for the propagating mode thus stands at a; = Q with the width F in 
the spectrum. The Vt represents the frequency of the collective oscillation. For 
the sound wave, for example, Vt = v^q where is the sound velocity. We note 
that because of the symmetry of the spectral function requested by general 
arguments, a propagating mode with the frequency VL is always accompanied 
by a propagating mode with the frequency — Jl, which propagates in the oppo- 
site direction in space. The emergence of a pair of propagating modes means 
that two slow variables must be involved. For the "particle mode," these two 
modes standing aX u = ±17 in the spectrum of course correspond to the par- 
ticle and anti-particle modes. As we approach the critical point, the width F 
gets small as a consequence of the critical slowing down. At the same time, 
the frequency Q. also goes to zero. This phenomenon is called the softening of 
a propagating mode. For the sound wave, the softening implies the decrease 
of the sound velocity. Noting that the sound velocity measures the stiffness of 
the system explains the meaning of 'softening.' 

It may be helpful to express the above statements in terms of the pole position 
in the complex uj plane. The diffusive mode is a mode that is located on the 
imaginary axis. Near the critical point, the pole moves on the imaginary axis 
toward the origin. The propagating mode is given by a pole located at an 
arbitrary point in the plane except on the imaginary axis, which shifts toward 
the origin as the critical point is approached. 

Thus, both the width F of the slow (either diffusive or propagating) mode 
and the frequency Vl of the propagating mode become small near the critical 
point, representing the critical slowing down and the softening, respectively. 
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Concerning how F and Q go to zero, Hohenberg and Halperin have proposed 
the dynamic scahng hypothesis [19]: Let Tg be a typical inverse time scale for 
each slow mode. We can take the width for Tg for the diffusive mode. For the 
propagating mode, the frequency or width plays the role of F^. The dynamic 
scahng hypothesis states that F, scales as a power law of the correlation length 

r,(0 = r7(90, (2.1) 

where / is some function. The index z is called the dynamic critical exponent 
and characterizes how the motion of the critical fluctuation gets slow or how 
the propagating mode is softened. The exponent z may differ from mode to 
mode at one critical point. The original dynamic scaling hypothesis insists 
that the exponents should be identical for all the modes in the system. This 
seems to be a reasonable assumption because the couplings among the modes 
are so strong that all the modes are in the same motion. However it is known 
that, in many systems a single value of z is not always shared among all the 
slow modes. 

As is already described, the non-equilibrium processes near a critical point are 
dominated by the long wavelength fluctuations of the slow variables. It is rea- 
sonable to expect that the idea of universality, which means that the dynamics 
is irrelevant to the microscopic interactions in the system, is still meaningful 
for the dynamic case just as for the static case. As is well known, the static 
universality class is determined solely by the symmetry of the system. The 
dynamic universality class, on the other, is dependent on what kinds of slow 
variables are contained, besides on the symmetry. Hohenberg and Halperin 
have defined their dynamic universality class [14]. The rules for classification 
are; (i) whether or not the order parameter is conserved, (ii) what kinds of 
conserved quantities are possessed by the system. According to the rules, the 
ferromagnet and the antiferromagnet, for instance, belong to distinctive dy- 
namic universahty classes because the slow variables in the two systems are 
quite different, although their static universality classes are the same. 

2.2 Dynamic universality class of the chiral phase transition 

Now let us consider the chiral phase transition which is of our major interest. 
The static universality class of the two flavor chiral phase transition is, as 
is well known, the same as that of the 0(4) Heisenberg ferro- and antiferro- 
magnet [1,7]. In Ref. [7], Rajagopal and Wilczek have discussed the dynamic 
universality class of the chiral phase transition. They argued that the chiral 
phase transition belongs to the same dynamic universality class as the Heisen- 
berg antiferromagnet, by comparing the slow variables of the both systems. 
The slow variables of the Heisenberg antiferromagnet are composed of the 
staggered magnetization N which is the order parameter and not conserved. 
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the magnetization M which is the generator of rotation, and the energy E. In 
the chiral phase transition, on the other hand, wc have the four component 
meson field 0^ which plays the role of the non-conserved order parameter, the 
six generators of the chiral rotation Q„, the energy E and the momentum P*. 
We note that the momentum is lacking in the Heisenberg antiferromagnet in 
which the spin is sticked on the lattice sites. The difference, however, is not 
so crucial in the following discussions. We have only to consider the itiner- 
ant antiferromagnet rather than the Heisenberg antiferromagnet for complete 
correspondence. Thus, according to Hohenberg and Halperin, the two systems 
belong to the common dynamic universality class and we could expect the 
same dynamic critical behavior in both systems. Since the antiferromagnet 
had already been well studied [20], the analysis went on in the same way. 
Rajagopal and Wilczek employed the kinetic equation provided by the time 
dependent Ginzburg Landau theory. By applying the renormalization group 
method to the kinetic equation, the critical exponent z — d/2 — 3/2 was 
obtained, where d is the spatial dimensionality. 

It is at this point where we should find the crucial difference between the 
antiferromagnet and the chiral phase transition by thinking of the slow modes 
rather than the slow variables in each system. Consider the disordered phase 
of the two systems. For the antiferromagnet, the slow mode associated with 
the order parameter is known to be diffusive. On the other, the slow mode 
for the chiral order parameter is apparently the meson mode or a "particle 
mode" [6,8,10], which is a propagating mode. This difference should be taken 
seriously. 

As is already mentioned, in general, wc need two slow variables in order to 
describe a propagating mode. This is a reflection of the fact that a propagating 
mode always accompanies its counter partner in the spectral function. The 
counter partner for a "particle mode" is just an anti-particle mode. What are 
the two slow variables for the description of the meson mode? The answer 
is that the two slow variables for the meson mode are the meson field 
which is the order parameter, and the canonical momentum tTq conjugate to 
This statement is just what we would like to propose and emphasize in 
this paper. Note that the canonical momentum is neither an order parameter 
nor a conserved quantity. Since the slow modes in the antiferromagnet and 
in the chiral phase transition are different, and accordingly, since the slow 
variables in the two systems are different, it is not possible any longer to 
identify the dynamic universality class of the chiral phase transition with that 
of the antiferromagnet. In fact, as we will calculate later, the dynamic critical 
exponent for the chiral phase transition will be found to be 2; = !— 77/2 = 
0.98 within our framework, which is apparently not the same a,s z = d/2 = 
1.5 of the antiferromagnet. The difference of the dynamic critical exponents 
explicitly shows the different dynamic behavior in the two systems. 
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We have found that the dynamic universahty class of the chiral phase tran- 
sition is not identical with that of the antiferromagnet. This means that the 
classification method of the dynamic universality class due to Hohenberg and 
Halperin is not apphcable at least to the chiral phase transition. What is the 
problem? The classification method of Hohenberg and Halperin is essentially 
to compare the slow variables in the system. The usual prescription to de- 
termine the slow variable, which has been accepted without any doubt, is to 
gather the order parameter and the conserved quantities of the system. As 
we have noted, however, the canonical momentum tt^ which is necessary for 
a proper description of the meson mode is neither the order parameter nor a 
conserved quantity. The chiral phase transiton presents us with the situation 
in which the order parameter and the conserved quantities do not exhaust 
all the slow variables necessary for a description of the slow modes. The pre- 
scription of just gleaning the order parameter and the conserved quantities 
for determination of the slow variables has thus broken down. Accordingly 
the classification method of Hohenberg and Halperin turns out insufficient 
because it involves just comparing the order parameter and the conserved 
quantities. If we consider all the slow modes or the slow variables necessary 
for their description, it can happen that the dynamic universality class given 
by Hohenberg and Halperin's method splits into finer classes. 

Now what we have to do is to reanalyze the chiral phase transition, the dy- 
namic universality class of which is not the same as that of the antiferromagnet 
nor that of any other critical point that has ever been considered. For that 
purpose, we will employ the mode coupling theory. In the analysis, we will find 
that a propagating mode for the meson mode appears appropriately from the 
meson field and the canonical momentum conjugate to it. In the next section, 
we will briefiy review the mode coupling theory. 



3 Mode coupling theory 



The mode coupling theory is a phenomenological theory which provides us 
with the kinetic equation for the critical fluctuations. It was developed and 
established by Kawasaki [15,16,17,18], and played an important role in the 
extensive studies of the dynamic critical phenomena in 1960-70's. Although the 
theory does not go beyond phenomenology, it gives us the useful quantitative 
predictions. 

In Subsec. 3.1, we will explain the general form of the kinetic equation in the 
mode coupling theory. In Subsec. 3.2, as a demonstration, we will give the 
kinetic equation for the Heisenberg antiferromagnet in the disordered phase. 
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3.1 Kinetic equation in the mode coupling theory 



Our kinetic equation is based on the Mori equation [21]. In general, it is nec- 
essary to separate the macroscopic degrees of freedom (slow motions) from 
the microscopic ones (rapid motions) in order to obtain a kinetic equation. 
This can be implemented by the projection operator method [22]. The Mori 
equation is derived by applying the projection operator method to the micro- 
scopic equation of motion, ie., the Heisenberg equation. The Mori equation 
deals with only the slow motions directly. In a sense, the Mori equation is just 
the Heisenberg equation after the rapid fluctuations are integrated out. Thus 
we can state that the kinetic equation in the mode coupling theory originates 
from the microscopic Heisenberg equation. 

The explicit form of the kinetic equation in the mode coupling theory is given 

by 



+^ E E ^rM (Xix?-.-)"' - ) + (3-1) 



where 



j If 



(3.2) 



eq 



2:(K-.<])«,4--'«-.4'vy • (3-3) 



The is the microscopic expression for the Fourier transformed fluctuation 
of the jth slow variable. The x| is the static susceptibility given by 

X|-<^^=(4'«0 = (4«^eq, (3.4) 

where the inner product {F, G) is defined by 

(F, G) = \ f d\l e^"Fe-^"G) (3.5) 

jj Jo ^ ' eq 

with the inverse temperature f3 = l/k^T and the Hamiltonian H. In the 
above expressions, (■ ■ • )cq denotes the statistical average in the equilibrium 
ensemble. In the following, we will omit the index 'eq' except for emphasizing 
it. 
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There are three terms in the right-hand side of Eq. (3.1), which are the hnear 
term, the nonhnear term and the noise term fi. The equation without the 
nonhnear term describes the free motion of the slow modes. The nonlinear 
term corresponds to the interactions among them. 

The hnear term consists of the frequency term and the dissipative term, which 
give the peak position and its width in the spectral function of the slow vari- 
ables, respectively. The frequency matrix ujji is given by the statistical average 
of the commutator of the slow variables in the equilibrium ensemble, as is 
shown in Eq. (3.2). For the classical system, the commutator is replaced by 
the Poisson brackets. By diagonalizing the frequency matrix, we obtain the 
normal coordinates out of the slow variables for the individual slow modes. 
The nonzero eigenvalue means the appearance of the propagating mode for 
the associated slow variables, while the zero eigenvalue corresponds to the 
diffusive mode. 

The L° in the dissipative term gives the Onsager (transport) coefficient. Within 
the present framework, L'^ is given by the correlation among only the short 
wavelength fluctuations, and thus does not show any anomalous behavior in 
itself. In general, however, the Onsager coefficient can be divergent at a criti- 
cal point because of the effect of the long wavelength fluctuation. In the mode 
coupling theory, actually, the renormalization coming from the nonlinear term 
can give rise to the divergence in the Onsager coefficient. In order to emphasize 
that L° is a bare or unrenormalized quantity, the suffix is put on L. 

The nonlinear term, which represents the interactions among the slow modes, 
is called the mode coupling term. It is known that the couplings among the 
fluctuations of the slow variables become very strong near a critical point [23]. 
Thus the mode coupling term is crucial for the critical dynamics. One of the 
effects of the mode coupling term is the possible divergence of the Onsager 
coefficient mentioned above. If the mode coupling term is dropped, the theory 
reduces to the van Hove theory [24] . 

We note that the coefficients in the kinetic equation are all given by the statis- 
tical average in the equilibrium ensemble and determined by the static prop- 
erties of the system. The dynamic nature is contained in the kinetic equation 
itself. Thus in the mode coupling theory, we input the static information into 
the kinetic equation in order to obtain the dynamic information. In particular, 
the criticality is furnished by the singularity of the static susceptibilities. The 
anomalous behavior at a critical point such as the critical slowing down and 
the softening is induced by the divergence of the static susceptibilities. One 
consequence is that the dynamic critical exponents are expressed in term of 
the static critical exponents in this description, as we will see later. 
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3.2 Kinetic equation for the Heisenberg antiferromagnet 

For illustration of the general formalism in the previous section, we here give 
the kinetic equation for the Heisenberg antiferromagnet in the disordered 
phase [15,16,17], which would be useful for the comparison with the chiral 
phase transition later. 

The Hamiltonian of the Heisenberg antiferromagnet is given by 

H ^ -J2 Jij^i • ^3 ' (3-6) 

where Si is the spin operator vector put on the ith lattice site which satisfies 

= ie-^^<5,,sl {a,(3,^ = x,y,z). (3.7) 

The Jij represents the exchange interaction and has a negative sign for the 
antiferromagnet . 

What we should do in the first place is to determine the slow variables. The 
slow variables are made up of the staggered magnetization (order parameter), 
the magnetization and the energy (conserved quantities). 

The Fourier transforms of the staggered magnetization (t|, the magnetization 
/i2 and the energy density are given by 

al^Y.r}'e-'^-'^{st-{st)). (3-8) 

i 

= Ee"''"'"^(«r-(0), (3.9) 

i 

H = E E i^i ■ - i^i ■ ^j)) ■ (3-10) 

i j 

where r)^ has an opposite sign for the nearest neighbors. 

The frequency matrix ujg- can be calculated by the commutation relation be- 
tween the spin variables. We find that the matrix elements are all zero. This 
indicates that all the slow modes are diffusive. Especially note that the slow 
mode for the staggered magnetization is diffusive as it should be. 

The dissipation terms can be found from the macroscopic hydrodynamics. 
Since the magnetization and the energy are conserved quantities, we have 

kBLl^k,,k/x,k-kX/x,j;, (3.11) 

^Bi^°,-,e.VXe.— ^'A7Q , (3.12) 



a P 
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which define the 'bare' Onsager coefficients. In particular, A° is the 'bare' 
thermal conductivity. The is the /c-dependent heat capacity per spin. For 
the staggered magnetization which is not a conserved quantity, we have 

^B^°,>,VX.,— i^°/x.,- ■ (3.13) 
The susceptibilities x's are defined by 

U^k^^-k)^ (3-16) 



where N is the total number of the spin site in the system. The coefficients 
of the nonlinear term, ^j-im, can be computed with the spin commutation 
relation. 

Finally we can obtain the kinetic equation. 



^ E f ^ - % X n,-^. + , (3-17) 



+ ^E(^-^)%x^^, + /e, (3.18) 

e.-= V^e,-+/|. (3.19) 



dt " " Cg 



Note that the energy mode decouples from the magnetization and staggered 
magnetization modes. 



4 Slow modes near the chiral phcise transition in the 0(2) sigma 
model 



In this section, we apply the mode coupling theory to the linear sigma model 
for investigation of the chiral phase transition. We examine what kind of slow 
modes, diffusive or propagating modes, appears near the critical point. For 
that purpose, it is sufficient to consider the frequency matrix uj in the kinetic 
equation. If the matrix element gives a finite value, it means that there appears 
a propagating mode in the associated slow variable. The zero matrix element 
indicates a diffusive mode. For simplicity, we consider the 0(2) linear sigma 
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model. The Lagrangian is 



02^)', (4.1) 



where 0a with a = 1, 2 is the meson field, and // and A are the bare mass and 
coupUng constant respectively. For dcfinitcncss, we assume that the 0i takes 
a finite vacuum expectation value = in the ordered phase. Thus the 
fluctuations of 4>i and 4>2 correspond to the a and tt meson modes, respectively. 

As we argued in Subsec. 2.2, the slow variables for the chiral phase transition 
are composed of the following; the meson field <i){x,t) (order parameter), the 
canonical momentum 7r{x,t), the chiral charge density Q{x,t), the energy 
density E{x, t) and the momentum density P*(x, t). The last three of these slow 
variables are the conserved quantities. Note that there is only one chiral charge, 
which is contrasted with the six charges in the 0(4) case. The microscopic 
expressions for Q{x, t), E{x, t), P'(x, t) are given by 



Q{x, t) = (7ri02 - 7r20i) (f , t), 
E{x,t)= I7ra2 + l(v0a)% 

P\x,t) ^ - (naV'cPa) {X,t). 



A 

+ 4 



(f, t), 



(4.2) 
(4.3) 
(4.4) 



Here, as usual, we expand (f)a{x, t) and 7Ta{x, t) into the normal modes in terms 
of the creation and annihilation operators. 

k y ak 
^ k 

where V is the volume of the system. The dispersion relation uo^^ is generally 
complicated as a consequence of the interaction. For the free system (A = 0), 

the Klein-Gordon equation gives ujj: = \fW+~jj?. We will determine at the 
final stage of consideration by imposing a consistency condition. 

The usual quantization conditions 

[(f)a{x, t) - {(f)a) , TTbiy, t)] = iSabS{x - y) , (4.7) 
[(t>a{x, t) - {(t>a) , 0(,(y, t) - (06)] = [7r„(f , t), 7r;.(y, t)\ = 0, (4.8) 
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or 



«afc(^)'«L-,W] =Sab6rkn (4-9) 

%kit)^ %k'it)] = Hkit)^ aK{t)] = 0, (4.10) 
provide us with the algebra to calculate the frequency matrix. 

The Fourier components of the fluctuation of the slow variables are given by 



Qq{t) = ^ab 



(4.11) 
(4.12) 
(4.13) 



+ 



A 



— (0a^) ] (^, t) — (the same expression)^^ 



(4.14) 
(4.15) 



where is the total antisymmetric tensor with ei2 = 1. Equivalently we have 



(t>aq{t) 



O'aq + ^a-q 



T^aq{t) =(-i) (aa^{t) - 



(4.16) 

(4.17) 
(4.18) 



1 1 -1 

E^V') ~2~Jy ^ ~^y^ak'^aq-k y^ak 
k 

+^^E[-^-(5"'-^)V 



"■Q-k ak- 



1 



aq—k 



A 1 ^ 1 
+ 7771 E J, 



4\/i 4 9.fci+-+fc4 . 

ki—k4 V 



^ak[^ak2^bk3^bki 



— (the same expression ) 5g-^o , (4.19) 



o-q—k ak—qj 



(4.20) 



The frequency matrix element is given by the commutator of each slow variable 
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in the mode coupling theory. By noting that the commutator between (pa and 
TTa is not zero, we can see a propagating mode appearing from these two 
slow variables. This propagating mode should and can be identified with the 
meson mode, which will be discussed later. Note that if we do not include 
the canonical momentum into the group of the slow variables just following 
the usual prescription for determination of the slow variable, this propagating 
mode would not appear and the order parameter fluctuation would undesirably 
show diffusive behavior. The canonical momentum plays an essential role in 
the chiral phase transition. 

The inclusion of vr^ is found to be reasonable if we consider the microscopic 
equation of motion, that is, the Heisenberg equation. The Heisenberg equation 
for the meson field consists of the two equations for the variables of and tt 
which are conjugate to each other; 

i^^cPix,t) = [Hx,t),H], (4.21) 
i^^n{x,t)^[7T{x,t),H], (4.22) 

where H is the Hamiltonian of the system. From these two equations, we 
obtain the Klein-Gordon equation for the meson dynamics. ^ Since our kinetic 
equation in the mode coupling theory is based on the Mori equation which 
is derived from the Heisenberg equation, it is natural that we need tTq as the 
degree of freedom for a description of the meson dynamics. 

Turning to the case of the antiferromagnet, we realize that the Heisenberg 
equation for the spin variables is 

ds" 

i^^[sf,H] {a^x,y,z) (4.23) 

with H being given by Eq. (3.6). Thus the staggered magnetization a = 
{a^,ay,a^) corresponding to the spin variable s = (s^,s^, s^) should be the 
degree of freedom for the dynamics of the antiferromagnet. In fact, from the 
variable of the staggered magnetization arises a diffusive mode as it should 
be, as we have already seen in Sec. 3. 

The canonical momentum as a slow variable is itself nothing new in the chiral 
phase transition. We know at least two examples of critical points in which 
the canonical momentum enters into the member of the slow variable. 



^ We note that the Klein-Gordon equation is of second order with respect to the 
time derivative while the Heisenberg equation is of first order. In order to reproduce 
the Klein-Gordon equation from the Heisenberg equation, it is inevitable to set up 
the simultaneous equations for two variables. 
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One is the A transition in the superfluid. The order parameter of the tran- 
sition is given by the wavefunction of the Bose condensate $. Because the 
wavefunction is a complex variable, the order parameter consists of the two 
components, i.e., $ and The wavefunction $ obeys the non-relativistic 
Schodinger equation, and as is well known, the canonical momentum conju- 
gate to $ is given by From the variables of $ and (and the entropy 
density), there appears the second sound mode which is a propagating mode 
(and one diffusive mode) in the superfluid transition. 

The other example is associated with the spin wave mode in the ferromagnet. 
In the ordered phase, supposing that the spontaneous magnetization points in 
the z direction, the transverse fluctuations of the magnetization, and My, 
make up the spin wave. We can regard the two variables of M± = ± iMy 
as conjugate to each other. 

Note that in the above two examples, the conjugate variable can be included 
into the member of the slow variable just by collecting the order parameter and 
the conserved quantity. In the chiral phase transition, however, the canonical 
momentum cannot be the slow variable by the prescription. This is nothing 
but the new feature of the chiral phase transition that we have not encountered 
in the systems so far considered in Ref. [14]. 

Alternatively to 0a and TTa, we can take 



for the two variables. These are obtained by dividing the meson field (paq-, and 
are just a liner transformation of 0a and TTa. In the following discussions, we 
will employ aa and a\ instead of 0a and tTq for the two degrees of freedom. 
This is because aa and give rise to the diagonal frequency matrix as we will 
see later and are more convenient. 

Thus for the two degrees of freedom to describe the meson mode, we will adopt 



1 




and 




(4.24) 



<l>aq{t) 



1 



(aaq- - (aaq-)) , 



(4.25) 





(4.26) 



Now we calculate the frequency matrix ujji given by 




(4.27) 



for the slow variables of 




(4.28) 
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where j on Aj specifies each slow variable. The static susceptibilities for (f), 
and (f>l_g' are given by 

Xaq = 



aq 



aqi yVaq j 



(4.29) 

(4.30) 



where the last equality holds because the susceptibility is a function of \q\^ . 
The other susceptibilities are defined similarly by 



XQq = {Qq, QI) , 

Xeq = [E^, -Si) , 
Xpq = (^Pqi 1 

where Xpq does not depend on i. 



(4.31) 
(4.32) 
(4.33) 



Firstly, let us consider the disordered phase, in which 0i and 02 are degener- 
ated. The frequency matrix is calculated to be 



where 



















































kYiT 




































'i^2qX2q 









































Xpq 














- 


—a' 

'Keq 





D 


= k^T 


k 






) Xak 





(4.34) 



(4.35) 



and the matrix elements are placed in the order of Eq. (4.28). We note that 
there appear propagating modes for 0i and 02, which should be identified with 
the a and tt meson modes. The kinetic equation without the dissipative and 
nonlinear terms is 



(t>aq = -i 



far 

'^^aqXaq 



(Paq, 



dt 



1^ (t'l-q 



'^^aqXaq 



(4.36) 



Here we impose a consistency condition that the frequency coincide with the 
dispersion ujaq, that is, 

^^aq/vaq 
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This is justified because in the mode couphng theory, the kinetic equation 
without the dissipation and mode couphng term should reproduce the Heisen- 
berg equation [16]. In fact, the equations given by Eq. (4.36) with Eq. (4.37) 
substituted are combined into 

^ (4.38) 

— (-i)a;„q-.(0„g--0^_-) = -a;„/(0„q.+ 0i_-), 

which give the Heiscnbcrg equation or the Klein-Gordon equation with the 
dispersion relation uja^. The condition (4.37) determines the dispersion relation 

as 



<^ag — 



In order to find Xaq-i consider the susceptibility for the field (f){x) itself. In the 
Ornstein-Zernike approximation, the susceptibility is written as 

(0(f), m') = ^ E w^/'''~''^ (^-^0) 

where m is the static screening mass. On the other hand, 

i^^m') -[iv^ (^^- + ^-.-) ^ ^ (4 + ^-.-') ^""'•') 

' ^2x^,e''<'-'^. (4.41) 



V 

k 



By comparing the two expressions, we find 



^2 _|_ 

Thus we obtain the dispersion relation 



Uaq = Z— 2^0^+771^. (4.43) 

We should note that the static screening mass m plays the role of the pole 
mass in this framework. As we approach the critical point, m goes to zero, 
which describes the softening of the meson modes. The way of the softening, 
or the way how m gets small, can be determined by the static scaling law, 
i.e., m ~ The renormalization constant Z yields the anomalous dimen- 
sion as Z ~ which leads to the dynamic critical exponent 1 — 7^/2 as is 
already advertised. Although the present argument does not take account of 
the nonlinear or mode coupling terms which are important near the critical 
point, the result does not alter even when they are properly taken care of. The 
full analysis will be given in Sec. 6. 
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dt ^ 


d pu 
dt^ 








= 0, 


dt 


by 




g2 J 





We can see that the slow mode for the chiral charge is diffusive. Moreover, the 
energy and the longitudinal component of the momentum couple to give rise 
to a propagating mode, while the transverse components give diffusive modes. 
The kinetic equations for the energy and the momentum are written as 

(4.44) 
(4.45) 
(4.46) 

where we defined and P?* by 

Pi = [5^^ - ^) PI + ^P| ^ Pf + Pf . (4.47) 

We note that here is the difference from the Heisenberg antiferromagnet in 
which the momentum is absent. In the Heisenberg antiferromagnet, the slow 
mode for the energy density is diffusive, while the chiral system gives the 
propagating mode for the energy because of the presence of the momentum 
density. What is this propagating mode? It would be essentially a sound wave, 
but not exactly the same as the usual first sound of the normal fiuid. For 
comparison, consider the slow modes in the normal fluid [17,25,26], in which 
the slow variables are the particle number density A^^-, the energy density E^, 
and the velocity u^^. The and may be transformed to the entropy density 
and the pressure through 

S^^^{E^-mhNi), (4.48) 

TTi 1 

Pq = ^^^^ - ^% (4-49) 

where p and h are the equilibrium values of the density and the enthalpy 
per imit mass, m is the mass of a molecule, and Xs(g) and as{q) are the 
g-dependent adiabatic compressibility and adiabatic thermal expansion coeffi- 
cients. The kinetic equation shows that the modes for and are diffusive 
while and combine to give the first sound mode. See Table 1. Thus our 
propagating mode is to be compared with the first sound in the normal fiuid. 
If the first sound is called the pressure wave, our propagating mode may be 
called the energy wave mode. If it is allowed to identify the chiral charge in 
the chiral system with the number density in the normal fiuid, we can say that 
the basis that diagonalizes the equation in each system is different. The basis 
of the one system can be obtained by the linear transformation from that of 
the other system. 

Now we move to the ordered phase. The frequency matrix in the ordered phase 
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is found to be 



UJjl 



V 

where 

















-B 

X.eq 








kBT 











B* 

Xeq* 
















A 












XQq 











kBT 


-A* 








2'^2fX2q 


XQq 








A* 


-A 

X2q- 











-B* 

Xlq 


B 

Xlq 














Xpq 























(4.50) 
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«10 + «lo 



(4.51) 



B =kT,T 



, 1 1 



-q,fel---fe4 



\kx^\k-}^hkz^hk^ 

t 



X 



'lfc2 



\-k2 



"bks 



a 



b-ks 



)( 



-k4 



(4.52) 



The 01 mode, i.e., the sigma meson mode can be treated in the same way 
as in the disordered phase. The 02 mode now couples to the chiral charge Q. 
We concentrate on the 3x3 matrix of |02g, Qq\- It can be diagonahzed 
readily to give the eigenvalues 



2\A\ 



\ \'2^2qX2q) X2qXQq 



(4.53) 



The zero eigenvalue gives a diffusive mode. The resulting propagating mode 
should be regarded as the 02 mode modified by the interaction or the coupling 
with the chiral charge. Thus we impose the consistency condition that the 



Chiral system 


Normal fluid 


Q 


Diffusive mode 




hi 


S 


■ • • Diffusive mode 


^1 


> • • • Propagating mode 






«\ 


> • • • First sound 




Diffusive mode 








• • • Diffusive mode 


6 1 



Comparison of the chiral system with the normal fluid. 
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eigenvalue be identical to ijj2q- Moreover since the 02 niode must be the Nambu- 
Goldstone mode, we should have a;2g = uq, where u is the pion velocity [8]. 
Thus we have 



\ \'2ui2qX2qJ X2qXQq 



uq. 



Furthermore if we assume 



2X2? 



then in the approximation of |g| -C 1, we find 

1 4M|2 1 



+ 



u^Z^ kBTxQq=o 



u 



(4.54) 
(4.55) 

(4.56) 



from which the renormalization constant can be determined if u is given. 
This treatment for the pion mode is essentially the same as that for the second 
sound mode in the superfluid transition [16]. 

We can see that in the ordered phase, the modes of the energy and momentum 
couple to the sigma meson mode through the nondiagonal elements. 

For summary, let us compare the ordered phase with the superfluid [16,26,27]. 
See Table 2. In the superfluid, the order parameter ($, $^") and the entropy 
S constitute the second sound wave mode and the diffusive mode. If we call 
the second sound the entropy wave, the tt mode may be viewed as the chiral 
charge wave in the analogy with the superfluid. 

Finally, we consider the hydrodynamic modes of the chiral system away from 
the critical point. In general, the conserved quantitiy and the transverse com- 
ponent of the order parameter in the ordered phase are slow variables even 
away from the critical point [26], the latter of which is slow because of a dynam- 
ical reason rather than a kinematical one. The fluctuations of those variables 
make up the slow (or hydrodynamic) modes and appear in the spectral func- 
tions as narrow peaks though their widths are finite for the finite wavelength 



Chiral system in the ordered phase 



Superfluid 



4>i 

02 

Q 

E 

L 



P 



a mode 
TT mode (NG mode) 
-I- Diffusive mode 

Energy wave 
Diffusive mode 




(Second sound (NG mode) 
-I- Diffusive mode 



P 



First sound 



Diffusive mode 



Table 2 

Comparison of the ordered phase of the chiral system with the superfluid. 
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fluctuations. In the chiral system, the conserved quantities, Q, E and P*, are 
the slow variables, so that the slow modes found in the above consideration 
will be responsible for the hydrodynamics of the system. In the disordered 
phase, since the meson fields (j)i and 02 are not conserved quantities, the fiuc- 
tuations of 0i and 02 do not make hydrodynamic mode and decay (or relax) 
in the microscopic time scale. Thus the a and tt meson modes will disappear 
away from the critical point [6]. In the ordered phase, 02 being the transverse 
fluctuation of the order parameter is the dynamically slow variable. Thus the 
TT meson mode appears as a hydrodynamic mode. The a meson mode, on the 
other hand, is not a hydrodynamic mode and decays in the microscopic time 
scale. This is, of course, consistent with the fact that we can not observe the 
sigma meson as a narrow peak in the vacuum. 



5 Kinetic equation for the chiral phase transition above Tc in the 
0(4) linear sigma model 

In this section, we give the kinetic equation for the 0(4) linear sigma model. 
For simplicity, we consider the disordered phase. The emphasis will be put on 
the effect of the propagating mode of the order parameter, which is replaced 
by the diffusive mode the itinerant antiferromagnet. 

The Lagrangian of the 0(4) linear sigma model is given by 

C^\{d,ct>S - \l^'<t>a' -\{<l>.')\ (5.1) 

where a runs from to 3 and 0^, denotes the meson fields, 0q, = (0o, 0i, 02, ^s)- 
For the pion field, we may use 0a = (0i, 02,03) with a running from 1 to 3. 
The conserved currents for the chiral SUl(2) ® SUr(2) transformation are the 
vector and axial vector currents given by 

= eaM''<Pc, (5.2) 

= 0„a>o - 0o5^0a, (5.3) 

where Cabc is the total antisymmetric tensor with £123 = 1. The chiral charges 
are then 

Q^{t) = J d'xV;^{x,t) = J d^xeabM<j)c, (5.4) 
Qtit) = J d'xAlix, t)^ j d?x {(PaMo - <Podt<pa) , (5.5) 

which may be accommodated in the antisymmetric tensor Qapit), i.e., 

Qapit) = / d^xeapjs(f>jdt(l)s, (5.6) 
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where eap-ys is the total antisymmetric tensor with 60123 = 1- The and 
are related to Q^p as 



Qoa — 1 Qab — ~^abcQ 



(5.7) 



The meson field 0^ and its conjugate momentum tTq are expanded in terms of 
the creation and annihilation operators, 



MS, ^) - 4^ E ^ i^akit) + a^_-(i)) (5.8) 

k V 



We note that the dispersion uoj: does not depend on a because all the a and 
TT mesons are degenerated in the disordered phase. We also note that the six 
charges, and , are degenerated in this phase. 

The quantization conditions 



[(f>a{x, t), (l)/3(y, t)] = [iTaix, t), 7r^(y, t)] = 0, 



(5.10) 
(5.11) 



or 



ait^ ^' (5k' 



= 



(5.12) 

(5.13) 



provide us with the rule to calculate the coefficients in the kinetic equation. 



The slow variables that consist of the Fourier components of the order parame- 
ter and its conjugate, the chiral charges, the energy density and the momentum 
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density are defined by 
1 



1 

Qal3q{t) — ea/375 



C^aq (^) ) 



V 

A 

+ - 



d^xe-'^-^(l)^Trsix, t), 

^ > (if, t) — (the same expression) 



(5.14) 
(5.15) 
(5.16) 



(5.17) 
(5.18) 



In the same way as Eq. (5.7), we can obtain Q^^ and Q^^ from Qapq by 

Qoaq — Q^q^ Q abq — — ^abcQ^- (5.19) 

The Qapq, Eg' and P| may be written in terms of the creation and annihilation 
operators as 



Qal3q{t) —^al3"/5 
Eq{t) 



1 



1 1 

'2 



E- 
-1 



1 



(5.20) 



+ \^l^[-H^-k)+A 1-;^= 
+^4e-^-^ ■ 



ak—q 



^ q,feiH |-fe4 



kl-kA V fcl ^2 fcs fc4 



— (the same expression) 5^^q , 



(laC + a. 



1 



V^2 

k 



1 ^. 



q—k 



0-ki) 

(5.21) 
(5.22) 



We shall now derive the kinetic equation for each variable, the general form 
of which is given in Eqs. (3.1), (3.2) and (3.3). 

We note that the second term in the expression of VLj-^un involves the frequency 
uj''J' and the presence of the propagating mode leads to appearance of the 
second term. Moreover the first term in Vlj.^im is also influenced by the presence 
of the propagating mode. Thus we expect a substantial difference in the kinetic 
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equation between the chiral phase transition and the itinerant antiferromagnct 
because the propagating mode associated with the order parameter is absent 
in the latter. 



The susceptibihties are defined by 



XQc^Qo:pq= [Qal3q, {Qal3q)^) = XQq , (5.23) 



Xeq^ {Eq^, {E^)^)^k^T^C^, 



Xpipiq= [^qA^g = Xpq 



where is the g-dependent heat capacity. Note that the susceptibihties are 
not dependent on the indices a or i — x,y, z. 

To obtain cui' and flj-im, wc need the commutation relations for the variables, 
which will be listed in Appendix A. 

Now we write down the kinetic equation for each slow variable. 
• Order parameter — (p^q 



'^X4>q, 



■"aq 



q, k, q-k 



Qa'p'q-k + ^ 



q—k 



^ak'^q—k 



P a'f3' 



+ fi , (5.24) 



where the dispersion cOg- is given hy cOg' = Z 2 ^^q-^ + in the Ornstein-Zernike 
approximation. We have defined the 'bare' Onsager coefficient for and 

4>l-q by 

kBLl^,aqlXaq = /^B^° tq-,atg-/Xatg- = i^J/^X-^g- (5-25) 
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The coefficient of the nonhnear terms, Vj-im = ^j;im/ ixiXm) ^ , are given by 



y I O: 0- Qa'ii' 

.q, k, q-k 



= -kBT I 



—1 



2VV 




{4>aqi ^^i3j;Qa'l3'k-q) 



^Qq-k '^^g'X<t>qX(l,kXQq-k 



(5.26) 



\q, k, q-kj ° ^ 



eq—k 



2uj. 



'qX^qX(j,kXeq-k 



, (5.27) 



\q, k, q—kj 



= -A;bT { 



Jap 



2VV 




, g, fe, q-ic 



= -A;Rr 



—1 



2v^ 





l—q 



{4>aq, 



"9' V^^fc-' k-q 



^pq-k "^^qXtpqX^kXp^-k 



(5.28) 




'9/ XQq-k 



{4>aq, 4>i3-kQa'f3'k-^ 



2u;. 



'qX<t>qX^kXQe[-k 



(5.29) 



ag*5 ^P-k^k-q 



'^'^qX<t>qX(f,kXeq-k 



(5.30) 



Vg, fc, g— A;/ 



Jap 



2VVV OJq 



CUT. 



^k i 
OJq 



^pq-k '^'^^XHX^kXpq-k 



(5.31) 
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For illustration, we shall give an explicit derivation of V ^' ! -) in Eq. (5.27). 



\q, k, q—kj 



{X(j)kXeq-k 

-kBT 



q—k 



^f3k k-q 
'^^qXH 



(5.32) 



The first term in the curly brackets vanishes. Note that the third term appears 
because of the presence of the propagating mode (the meson mode). For the 
second term, we use the approximation that 



^k^ Bp 



^fc+f' Bk/_i 



(5.33) 



where |A;'|, |/| -C 1, which has a ground in the locality of the microscopic 
interaction [17]. With this approximation, we can manipulate the second term 
in Eq. (5.32); 



(second term) = (0^- \(f)cq, E^__ 



13k 



yK'^ak^'f'U 




(5.34) 



where Ti. is the Hamiltonian and we have used the Heisenberg equation iO = 
[0,n] and the relation i(i,5) = kBT{[A,B]). Thus we obtain Eq. (5.27). 
We note that the second term of Eq. (5.32) is proportional to the frequency 
of the meson mode so that it vanishes if there is no propagating mode. Thus 
we see that the coefficient of the coupling of the order parameter and the 
energy density remains even in the disordered phase owing to the presence 
of the propagating mode of the order parameter. This is contrasted with the 
antiferromagnet, in which the mode of the order parameter is diffusive. In fact, 
in the antiferromagnet, the energy mode decouples from the other modes in 
the disordered phase as we have seen in Subsec. 3.2 [16,17]. We also note that 
the second term in Eq. (5.27) involving the three point correlation has arisen 
due to the propagating mode. 



Order parameter conjugate — 
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Similarly, we obtain the kinetic equation for 



a—q 



2^^^ (q, '%!q-^ " ) ^f3kQa'f3'q-k + ^ ^ak^q-k 



13 a'l3' 



(5.35) 



with the coefficients 



V 



q, k, q-k 



-kuTl 



2VV 



+ 



'^?/ ^Qq-k '^'^qX4>qX4,kXQq-k 



(5.36) 



\q, k, q—kj 



t J)t 

^Ok k-q 

'^'■^qX4>qX(j)kXeq-k 



(5.37) 



\q, k, q—k 



UJt 



J Xpq-k '^'^qX<l>qX,pkXpq-k 



(5.38) 



\q, k, q-k 



1 _^ {^a-q^^f3-kQa'l3'k 



(5.39) 



A;rT 



+ 



«-(?' ^l3-k^k-q 



eq—k 



2u, 



'qX<l>qX(j)kXeq-k 



(5.40) 
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Chiral charges — Qapq 



Since all the six charges should behave in the same manner, we take specifically 
Qiq — Qoiq- Its kinetic equation is given by 



d_ 



XQq 



q, k, q—k 



+ VI 

W, k, q-k 

with the coefficients 

y /Qoij 2, y /Qoi_, 2t,_3t\ 

\q, k, q—k) \q, k, q—k J 



k^3q-k'^ \q-^ q-k ) ~^2-fc~^3fc-q 

A niA 



Qoi, 2t, 3\ ,t A. /"Qoi, 2t, 3t\ ,t x\ 



^Ik^q-k 



(5.42) 



= -A;Rr 



q—k 



'q-k 



1 



2yy 



q—k 



+ 

fc V ^-f-^ ^<t>q-k, 



y /Qoi^, 2, 3_t\ ^ y /Qoi, 2t, 3 
Vq, fc, q-k J \ 



(5.43) 



V 
V 



= -A;rT 



q, k, q-k 
q, k, q-k 



^q, k, q—k 



2VV 



+ 



q—k 



UJ 



q—k 



y j Qoi, Q2' Q3 

\Q, k, q-k 



^k Xcf>k 

i / 1 



1 



^q-k X^q-k 
1 \ 



^ V^Qfe ^Qq-k, 



(5.44) 

(5.45) 
(5.46) 



We have defined the 'bare' Onsager coefficient Lq for the chiral charge by 

^^LQ^pq,Q^pql ^Q^m = I^Lq/XQq, (5-47) 
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which involves (f because the chiral charge is a conserved quantity. 
• Energy — Eg 



ka 



ka 



2 \<1, k, q—kj a—k ak—q 

ka 

+ iEv(j|;^:j)i5(Kj + /|, (8.48) 

ki 

The p and p in the dissipation term are the proper energy density and the 
pressure in the equihbrium state respectively. The A*^ is the 'bare' heat conduc- 
tivity and Cq' is the g^dependent specific heat. The dissipation terms for the 
energy and the momentum are derived from the relativistic hydrodynamics. 
The derivation will be given in Appendix B. 



The frequency is given by 



ep'q 



-eqyxpq, (5.49) 



and the coefficients are computed as 



V h ^' ! = kBTQq' ^ "' (5.50) 

\q, k, q—kJ " ^ -y ^-y ^ ^ ^ ' 

X-pqA(t,k^(t)q-k 



^^^^^^-^^ \ \/V '2uJkX^kX^cl-k XpqX^kX^q--k I' 



V h ff,) = kBTOq' ^^'' \ (5.52) 

\q, k, q-kj " ^ v -V -V - 

^' ' ^ I Xek V XekXpq-k XpqXekXp^-k 

(5.53) 



where 



^^li:{-i + \k')x,k- (5-54) 
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• Momentum — Pi 



^ - io; E A° ( ^ ^'^^"l 



+ 



+ 



ka ka 

2 V?, k, q-k J ^a-k^ak-q 2 4^ W, fe, q-kJ V'k^q-k \^k^q-k/J 

ka kO. 



k 
kjl 



(5.55) 



The ?7° and C° in the dissipation terms are the 'bare' shear and bulk viscosities. 
The dissipation terms are derived in Appendix B. 

The frequency is given by 



i^pieq = -©?VXeq-, 



(5.56) 



and the coefficients are computed likewise as 



\q, k, q—kj 
1 



q—k 



2VV^V 



-{q - kf 



q—k 



UJ 



q — k 



1 



^q-k 



^4>q—k 



eq 



aq—k 



XeqX(l)kXrf)^-k 

(5.57) 



Vq, k, q—k/ 
1 




(5.58) 
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q-k \V '^k \^q-k 



(5.59) 



\q, k, q—kj 



W V XQk 



(5.60) 



1 (^?. -g-g-gg-.-) 



(5.61) 



-Qq' 



XeqXpkX. 



pq—k 



(5.62) 



We note that there are time derivatives in the dissipation terms for Eg' and P|. 
These time derivatives may be ehminated by means of the kinetic equations 
themselves. We can perform it simply by replacing and P| with the asso- 
ciated frequency terms. This is because the time derivative of E^ and P| and 
the frequency term are of second order with respect to the small quantities, 
that is, the wavenumber and the fluctuation of the slow variables, whereas the 
dissipation term and the nonlinear term are of higher order. Thus we have 



2\0 



1 - 



T ^ 
p + p k-BT"^ 



(5.63) 



^ — Eq^+X'f^—\1 

Xeq Xpq \ 



T e ^ 

p + p k^T"^ 



Pq) 



V 







p+p \ ^ 3 / p+p 



(5.64) 



where the dissipation term and the nonlinear term are omitted. 



32 



6 Dynamic critical exponents 



We have derived the kinetic equation for the chiral phase transition in the 
previous section. For the full analysis, the numerical calculation would be 
involved. As an application without the numerical work, we shall estimate 
the dynamic critical exponents from the kinetic equation in this section. 
The dynamic critical exponents specify the way how the typical inverse time 
scales of each slow mode, i.e., the widths and the frequencies are scaled. The 
fluctuation-dissipation theorem may allow us to imagine that the characteris- 
tic scaling law of the dynamic fluctuation of the slow variable originates from 
the characteristic behavior of the fluctuation in equilibrium states. Then it 
is reasonable to expect that there are some relations between the dynamic 
and static critical exponents. Actually, the dynamic exponents are given in 
terms of the static ones within the mode coupling theory as we will see below. 
We also discuss the anomaly of the Onsager coefficients, which are caused by 
renormalization from the nonlinear terms. 

The analysis goes exactly in the same way as in Ref. [17]. In Sec. 6.1, we 
prepare the general formalism for obtaining the dynamic critical exponents. 
In Sec. 6.2, we apply the formahsm to the chiral phase transition. 



6. 1 General formalism 



In order to obtain the dynamic critical exponents, we have to know the char- 
acteristic scaling behavior of the dynamic fluctuations of the slow variables. In 
the mode coupling theory, the scaling properties of the dynamic fluctuations 
are given by those of the static ones as a starting point. The amplitude of 
the characteristic static fluctuation of the slow variable Aj is measured by its 
static susceptibihty Xj — (l^jP)- We define the "reduced" slow variable Aj by 

= ^j/VXr (6.1) 

We also rescale the kinetic equation (3.1) to obtain the reduced kinetic equa- 
tion for Aj, 

I \ \/XiXlJ ^ Ira 

(6.2) 
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where 




The reduced variables have a finite amplitude of order unity, so that the infor- 
mation of the characteristic scaling laws of the fluctuations of the slow vari- 
ables is now pushed away into the coefficients of the reduced kinetic equation, 
which are given in Eqs. (6.3) and (6.4). 

We consider the time correlation function for the reduced variable defined by 

gjiit) ^ {Aj{t)M0y) , {t>0) (6.6) 

which we will call simply the propagator. The physical information is contained 
in those propagators. The dynamics of the reduced variables is subject to the 
kinetic equation (6.2). By solving it, we can calculate the propagator. 

In order to analyze the propagator, it is convenient to introduce the second 
quantized formalism, in which we can make the field theoretical techniques and 
considerations. In the second quantized formalism, the propagator is written 
as 

gji{t) = {0\aje-''''ai\Q). {t > 0) (6.7) 



Uj and oij are the operators that create and annihilate the mode j. We have 
defined the associated "vacuum" state |0), in which no j modes are present. 
The "Hamiltonian" H governs the dynamics of the modes and is given by 

H = Ho + n\ (6.8) 

7^0 = - E hi + i (XiXO"^ 4] (6-9) 
jl 

7Y' = - - E {P'*j;lm^l^^(^3 + ^r,lm^3^l(^rn) , (6.10) 
jlm 

where Tfo and 7i' are the free and interaction "Hamiltonians" and correspond 
to the linear and nonlinear terms in the kinetic equation respectively. We note 
that the free "Hamiltonian" is not Hermitian because of the dissipation term. 
In the following, we will use Ti^^ and Tig^-^ to denote the two kinds of terms in 
Eq. (6.9), which describe non-dissipative and dissipativc processes respectively. 
In the second quantized formalism, the expression of the propagator Eq. (6.7) 
can have the interpretation that the mode I is created at i = 0, which turns 
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into the mode j through the mixing and couphng with the other modes to be 
annihilated at a later time t. 

We note that in this formahsm, the reduced coefficients cD's and Q's, which 
contain the information of the characteristic scaling properties of the fluctu- 
ations, appear in the "Hamiltonian" Ti,. Thus the "Hamiltonian" determines 
the scaling law of the slow modes. 

Our present goal is to obtain the dynamic critical exponents or the char- 
acteristic time scales of each mode. For that purpose, it becomes necessary 
to consider each mode separately and derive the effective "Hamiltonians" for 
each mode. Wc can obtain the effective "Hamiltonian" by introducing the pro- 
jection operator that separates the subspace of the interesting variables from 
that of the other variables. We denote the set of the slow variables for which 
we want the effective "Hamiltonian" by {a} and the eliminated set by {b}. 
The "vacuum" state is then a product of the "vacuum" states for {a} and {b} 
denoted by |0a) and \0b), 

|0) = |0„)|05). (6.11) 

The projection operator Vb is defined by P^l ■ ■ • ) = |0;,) (Ob| • • • ), which projects 
arbitrary states of {b} into the ground state |0b). 

Consider the one-sided Fourier transform of the propagator for the slow vari- 
ables j and I that belong to the set {a}; 

/■oo . \ 

g,i{u) ^ / dte^'g,i{t) = i{0\a, -&i\0) (6.12) 

Jo uj — rt 

We can manipulate it into 

gji{u) = i(0„|aj , . ai\^a) (6.13) 

UJ - na[i^) 

in order to have the frequency dependent effective "Hamiltonian" Ti,a{uj), 
which acts on the variables {a} with the interactions with the variables {b} 
included and is given by 

= Hoa + {Ob\n'\Qb) + {Ob\n' Vb) nr' (i - n) ^'lo;,). (6.i4) 

Now what we have to do for investigation of the slow modes is to know the 

scaling properties of the coefficients cD's and Q^s. We see that it is necessary 
and sufficient to have the scaling behaviors of the commutators of the reduced 
variables, i.e., [A, A], which we will examine for the chiral phase transition in 
the next subsection. 
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6.2 Investigation of the chiral phase transition 



Firstly, we consider the scaling behavior of the amphtude of the slow variable 
Aj itself. For the chiral phase transition, Aj consists of {0^^, 4>^^_j^, Qk, E^, P|}. 
These variables are the Fourier transforms of the corresponding density, which 
we will denote as Aj{x). The Aj{x) exhibits characteristic anomaly at the 
critical point and is scaled in terms of the correlation length ^. To make it 
explicit, we write Aj {x, ^, V) for Aj (x) where we have also included the possible 
dependence on the volume V. Applying the block spin transformation with 
the block size L on Aj{x, ^,V),we have the scaled amplitude of the fluctuation 
as 

Aj{x,^,V) = L-^^A.^x/mmV/L"), (6.15) 

where Xj is the exponent to be determined. For the Fourier component which 
is defined by 

A{k) = ^J d'^xe-"'-'^A{x), (6.16) 

we then have 

Aj{k, e, V) = L'^-^^Aj{Li i/L, VIL\ (6.17) 
The exponent Xj can be determined by the (static) susceptibility xf 

L''-'^^{\A,{Q,^mviL')\') = (i^,(o,e,^)r) = x,,g=o— e^. (e-is) 

where 7^ and v are the static critical exponents. When choosing L = ^, we 
find immediately 

For the chiral phase transition, we have 

1^ = 7, 7q = 0, 7e = a, 7p = (6.20) 

leading to 

A = ^, (6.21) 

(6.22) 

, 1 / N 

- =ci--, 6.23 

V ) V 

(6.24) 

where we have used the scaling relations 

q; + 2/3 + 7 = 2. (6.25) 



Xj^-[d-^]. (6.19) 























d 




xq = 


2' 








'd 








d 




Xp — 


2' 
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Wc note that the static critical exponents a,/3,7, v are the same as those of 
the ferro- and antiferromagnet because the chiral phase transition is in the 
same static universahty class as these systems. 

Thus the amplitudes of fluctuation Aj have the characteristic scaling proper- 
ties as follows; 

Qap{k,C,V)=eQap{^k,l,V/^''), 

E{k,^,V)=^^E{^k,l,V/e), (6.26) 



We now discuss the scaling property of the commutators of Aj. If the exponent 
of [Aj, Ai] is denoted as Xji, namely [Aj, Ai] scales as in the same sense 

as Eqs. (6.26), then we find that [^j,^;] scales as ^^i-^^i-^l'^-^n ^ that is, 

(6.27) 

The computation of Xji for {0^^, 4>\_^') Qk^ E^, Pi} in the chiral phase transi- 
tion is performed in Appendix C. Using those results, we find 

[0,0t]^r-, [0,0] = [0t,0t] = o, 

[4(5]~ri [Q,^]~r^-^ [Q,p]~r^-\ (6.28) 

[p,p]^r'^-\ 



We are now ready to derive the scahng behavior of the coefficients in the 
kinetic equation or in the "Hamiltonian." The nonvanishing frequencies of the 
propagating modes are those of the meson mode and the energy wave mode. 
The scaling properties for the reduced frequencies are 



-kBT 
-k^T 



t 



For the scaling properties of the Onsager coefficients, we have 

tO 7-0 t-2 \0 ^0 aO t-2 



(6.29) 
(6.30) 

(6.31) 



because (j) is not conserved while Q, E, and P* are conserved. We also note 
that 

X</.~C^, XQ~Xp~C°, Xe^i'^- (6.32) 
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Thus we can find the scahng properties of the free "Hamiltonian" in the second 
quantized formahsm. The free "Hamiltonian" Hq is separated into several 
parts: 

rto = rto<l> + riQcl> + rlQQ + ^Oep^ + «0e + ^OpT + /topL, [p.66) 

the scaling properties of which are found to be 

'^oQ~C~^ 

<ep^-r^-^-^ ^o^e-r^-^, (6-34) 

'I/O t-2 'I/O , /:-2 

We next turn to flj-im- This is contained in the interaction "Hamiltonian" 7i' 
and determines its scaling property. We note that flj^im involves the commu- 
tators of Aj. Hence it would be useful to divide the interaction "Hamiltonian" 
as 

n' = % + n'Q + K + n'p (6.35) 

with 

'^e = ^[e,<I)] + ^'£;,Q] + ^[e,e] + ^[E,p]) (6.36) 
where the scaling properties of each term are 

(6.37) 

Wc are now in a position to discuss the dynamic behavior of each slow mode. 
We start with the dynamic behavior of the meson mode or the order param- 
eter fluctuation. The information of the dynamic behavior is contained in the 
reduced propagator 

9ct>ga(t) = (0|Q;^g-ae-'*^a^g-„|0), (6.38) 

where a^^a and a^^a are the annihilation and creation operators for the mode 
(paq- Prom Eqs. (6.13) and (6.14), we find that the Fourier transform of g^^a{i) 
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is written as 

g(l>qa{<^) = i(O0|Q;</,ga ^ / x Q:^gh|0^), (6.39) 

u — n(f,[u!) 

where 

7^0(0;) + + (OQep|H'|OQep) 

+ {OQep\n' ^ (1 - T^Qep) ^lOgep). (6.40) 

The vacuum states jO^) and \OQep) are defined with respect to the modes {0} 
and {Q, e,p} respectively and Vqep is the projection operator onto |OQep)- 

Here we need to use the general property of the interaction "Hamiltonian" [17]. 
The general argument tells us that the interaction "Hamiltonian" involving 
the commutator of a mode j has at least one creation or annihilation operator 
of j. This means that 

m%,i]\Oj)^0. (6.41) 

Moreover it is proved that the interaction "Hamiltonian" with flj-im, ^ e., H'j 
has at least one creation operator as well as one annihilation operator of j. 
Thus both n'jlOj) and {OjlUj vanish. 

With these in mind, we find 

(Ogepl^'IOQep) = {OQep\n[^JOQep) ~ T^- (6.42) 

As for the last term in Eq. (6.40), we note that at least one mode among 
{Q,e,p} must be excited in intermediate states because of the presence of 
1 — Vqcp- In order to make the last term as large as possible, it is desirable 
that the denominator is small. Because the "Hamiltonian" for the mode, 
~ '^[cj) 4>] ~ 1 is larger than that of the other modes, the denominator 

would become larger in the absence of the mode in the intermediate state. 
However this is not allowed because 

(OQep|H'|0<^) - (0^|H'|OQep) = 0. (6.43) 

In the end, it turns out that the largest contribution is found in the presence 
of the (f) mode in the intermediate state to be (?^|^ ^j)^/ (T^^ + Hj^ ^j) ~ ^-^l"^^ . 
Thus the two interaction "Hamiltonian" s in the effective "Hamiltonian" for 
the mode arc determined by the mode itself. This indicates that the 
mode almost decouples from the other modes. 

We have clarified the scahng property of the effective "Hamiltonian" for the 
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mode, 

+ (pQepI^' ^ _ ^-^ _ (1 - Vqep) T-C'lOQep) 

.^^-5^7+^-5+^-517+^-^. (6.44) 

The dynamics of the (p mode is dominated by TCq^ and the two interaction 
terms, all of which scale as ^"^1"^^ . This means that Qt^qa {t) has the following 
form, 

94>?a{t) = F(tT^q^^{0:QO (6-45) 
with the characteristic frequency of the form, 

r^,-a(0 = #/(90- (6-46) 
Thus we obtain the dynamic critical exponent z for the mode 

^*=^ = l-|. (6.47) 

When we use the values of the static critical exponents, rj = 0.03 [28], we find 
Z(j, = 0.98, which is contrasted with the value zn = d/2 = 3/2 of the staggered 
magnetization mode in the antiferromagnet. 

Recently Boyanovsky and dc Vega have calculated the dynamic critical ex- 
ponent Zfj) in a different approach [9]. They apply the renormalization group 
method directly to the microscopic theory, i.e., the field theory of the O(A^) 
linear sigma model and obtain 

For = 4 and e = 1, the above expression gives z^p = 1.04, which is compatible 
with our result rather than the value for the antiferromagnet. 

We note that the fiuctuation effect reduces the value to less than unity in our 
framework, while the opposite occurs in the calculation in Ref. [9]. The further 
close investigations of the kinetic equation, i.e., the appfication of the dynamic 
renormalization group method to the kinetic equation, may reveal the critical 
exponent larger than unity. Or it might be possible that the inclusion of higher 
order terms in the e expansion and the proper resummation of those terms in 
the approach in Ref. [9] turn its tendency to the opposite. In any case, it needs 
more involved investigation to settle whether the dynamic critical exponent is 
larger or less than unity. 

We note that if we employ the values in the large N limit for rj, that is, f] = 0, 
we have z^f, = 1. This is reasonable because the fiuctuation is suppressed in 
the large N limit. 
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We find that the dissipation term is smaller than the interaction terms. 
Thus the dissipation term is overwhelmed or renormalized by the interaction 
terms, which results in the same scaling pattern as the frequency term, that 
is, ~ As a consequence, the Onsager coefficient diverges at the critical 

point as 

L^r^i^. (6.49) 

What this means is that the width of the meson mode becomes narrower just 
in the same manner as that of the softening when we approach the critical 
point. This is not necessarily trivial. In fact, one of the motivations of the 
calculation in Ref. [9] is based on the suspicion that the meson mode may 
not be a good quasiparticle mode with a narrow width near the chiral phase 
transition. The calculation shows that the meson mode indeed becomes a good 
quasiparticle near the critical point, which agrees with our result. 

As a final comment on the mode, we give the comparison with the antiferro- 
magnet. In the antiferromagnet, the scaling property of the commutator of the 
staggered magnetization a is given by [a, a] ~ ^"'^^^ which leads to Ti.' ~ i^'^^'^- 
The dissipation term scales as ~ so that it is renormalized by 

Ti' to scale as ~ ^-'^Z^. This gives the critical exponent — d/2, which is also 
predicted by the renormalization group method [20]. 

Next we consider the chiral charge Q. Using Eqs. (6.13) and (6.14), we have 
the propagator of the Q mode as 

gQqanii^) = i(0Q|Q;Q^a/3 }, , x Q:Qga^ |0q), (6.50) 

where 

+ (^t>-p\'^' ^-^ilr^^^)n ~ ^^''^ ^'\^^ep). (6.51) 



The second term reduces to 



pep\H[Q^Q-\\Q^ep) ^ i 2. (6.52) 



The similar argument to the (j) mode dynamics applies to the third term of 
Eq. (6.51). There arises the most important contribution from the situation 
where the intermediate state has the Q mode as well as the e or p mode 
excitation but does not have the (f) mode. Because the Q mode "Hamilto- 
nian," Tijg Qj ~ is greater than those for the e and p modes, the Q 
mode dominates in the intermediate state. Therefore the third term scales 
as (^[(3,(5])^/^(q,(5] ~ C"*^^^- As in the (p mode, the two interaction "Hamilto- 
nian" s in the effective "Hamiltonian" Hq are determined only by the Q mode. 
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which means that the Q mode decouples from the other modes. Thus we find 

1 



— ^OQ + (0</>ep|'H'|0</,ep) + {Otf,ep\Ti-' _ /-, _ ^ 7^77 (1 ~ 'Pcpep) '^'\^(t>ep) 



~r'+r^+ri (6.53) 



which gives us the dynamic critical exponent for the Q mode, 

ZQ='^-. (6.54) 

We see that the dissipation term receives the renormalization from the inter- 
action terms and the Onsager coefficient diverges as 

Lq^e-i- (6.55) 

The behavior of the Q mode is quite the same as that of the magnetization fi 
in the antiferromagnet which is the counterpart of the Q mode in the chiral 
phase transition. The magnetization is renormalized by the interaction terms 
to have the critical exponent — (i/2, which is the same as that of the 
staggered magnetization. In this sense, the Q mode may be regarded as a 
vestige of the antiferromagnet. 

The third mode to be considered is the mode. The discussion goes in the 
similar way to before. The propagator for p^ is written as 

gp-^qii.^) = i(OpT|Q;pTq-i ^ . x ^pTg-^lOpT), (6.56) 

CtJ — rlpTyUJ) 

where 

^pT(a;) =HqpT + {^4>QepA^'\^4>Qep^) 

+ (O^QepLl'H' J- — ^ (1 - Vcj^Qep^) 'H'lO^Qepi^). (6.57) 

- (^-L - r<t>Qep^ ) ri 

The scaling property of the second term in T-CpT{u!) is given by 

{^<l>QepAT^'\^<t>Qep^) = {^<l>QepA'^'lp^ ,p'^]\^ct>Qep^) ^ (6.58) 

For the third term, we have the most important contribution coming from the 
process with the p"^ as well as the e or p^ modes but without the and Q 
modes in the intermediate state, which scales as ~ ^-d+a/u^ Thus we find 

HpT — Ti-QpT + (O^Qep^Ti'lO^Qep^) 

U — (1 — V^Qep^ j ri 

-r'+r^-^+r'^^. (6.59) 
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In the intermediate state in the third term, the e and modes couple to the 
mode. The most dominant term in the effective "Hamihonian" is, however, 
not the interaction terms but the free term. Thus the mode again closes its 
dynamics within itself. The mode scales as ~ without renormalization, 
leading to the critical exponent 

ZpT = 2. (6.60) 

Since the renormalization is not effective, the Onsager coefficient of the trans- 
verse momentum, i.e., the shear viscosity, does not diverge at the critical point. 
This may be compared with the transverse velocity mode in the normal fluid 
in which the shear viscosity shows the logarithmic divergence. The logarithmic 
divergence is a consequence of renormalization from the third term with the 
^"^ scaling. In the case of the normal fluid, there is a very slow mode, that 
is, the entropy mode which scales as Because the intermediate state is 
dominated by this entropy mode, the denominator of the third term becomes 
very small, which results in the larger third term than ours. 

Finally we consider the energy wave mode composed of the e and variables. 
The propagator is found to be 

9ep^qi{^) = i(OepL|aepLg-j — ;^Q;epLg-i |OepL) , (6.61) 

where 

+ {^ct>Qv'^\'X-' 7- ^ (l - V^QpA n'\O^QpT). (6.62) 

UJ - [l- V^QpT jH^ 

The first interaction term scales as 

(O0QpT|7i'|O^QpT) = {0^qpt\H[e^e] + 'H[e,p^] + '^[pL,pL]|O0QpT) ~ 

(6.63) 

The most important contribution of the second interaction term arises from 

the process in which the e and p^ modes as well as the p"^ mode appear but the 
and Q modes are absent in the intermediate state, which gives the scaling 
behavior of the last term ~ ^<^-(4-a)/!/_ ^jj-^^g 

Hepi^ = "^O^pL + Hq^ + HopL + {O^QpT\Ti'\O^QpT) 

+ (0«pT|?i' J- (l - V^QpTj n'\O^QpT) 

Lj- [1- v^QpT jn^ 

where the two terms in the parenthesis correspond to the terms associated 
with the thermal conductivity A° and the bulk viscosity C° in "^^l- We see 



43 



that the non-dissipation term is renormahzed by the interaction terms, but 
the dissipation terms are the dominant terms and not renormahzed, which 
give the critical exponent 



^epL = 2 + - ^ 1.74 (6.65) 



with a — —0.19 and u — 0.73 substituted [28]. Thus the timescale of the 
energy wave mode is determined by the e and modes themselves and is 
not affected by the other modes. We note that the bulk viscosity does not 
diverge in contrast with the normal fluid. In the normal fluid, the commutator 
[p, v}"], that corresponds to [E, P^] in our system, is very large. It enters into 
the non-dissipation term as well as into the latter of which thus renor- 
malizes the dissipation term. Because the slow entropy mode dominates the 
denominator, the second interaction term becomes even larger. Thus the bulk 
viscosity of the normal fluid shows the very strong divergence. In the chiral 
phase transition, [E, P^] is somehow very small and the intermediate state 
in the second interaction term is not occupied by such a slow mode as the 
entropy mode.Thus the divergence of A and C, does not occur. 

Prom the above discussions on the four modes, we notice that the dynamic 
critical exponents, i.e., the time scales of each mode are all different and all 
the modes decouple from each other. As mentioned in Sec. 2.1, the original 
dynamic scaling hypothesis proposed in Ref. [19] claims that all the expo- 
nents are identical. In this sense, the "original" dynamic scaling hypothesis is 
violated in the present system. 

Although the four modes decouple completely in our system, this should be 
regarded as an extreme case. In the antiferromagnet, the modes of the stag- 
gered magnetization and the magnetization couple to each other strongly. The 
critical exponents for the two modes are identical, though the energy mode 
has a different critical exponent. Thus the original dynamic scaling hypothesis 
holds partially in the antiferromagnet. 

We note that in the chiral system, the largest critical exponent or the largest 
time scale is provided by the mode. This means that the mode is the 
slowest among the four modes. Thus it determines the relaxation time of the 
whole system. The other modes in themselves are certainly of interest. In 
particular, the softening and the narrowing of the meson mode have the many 
physical implications. But when we consider the relaxation of the system, it is 
dominated by the decay of the p^ mode, and the other modes are not relevant. 
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7 Summciry 



We have discussed the dynamic aspect of the chiral phase transition. The 
classification method of the dynamic universahty class had been established 
by Hohenberg and Halperin. The method is based on the prescription of how 
to collect the slow variables in the system. The prevailing prescription tells 
us to take the order parameter and conserved quantities of the system for 
the slow variables. If we follow the prescription, we are incorrectly led to 
identify the dynamic universality class of the chiral phase transition with that 
of the antiferromagnct and to find the dynamic critical exponent of the order 
parameter fluctuation to be 2; = d/2 — 1.5, as Rajagopal and Wilczek had 
argued. 

We have then turned into the consideration of the slow mode rather than the 
slow variable in order to find the crucial difference between the chiral phase 
transition and the antiferromagnct. While the order parameter fluctuation of 
the antiferromagnct in the disordered phase is a diffusive mode, the meson 
field which is the order parameter of the chiral phase transition gives appar- 
ently a propagating mode. We have stated that in order to describe the meson 
mode appropriately, we must include the canonical momentum conjugate to 
the meson field into the member of the slow variables. Since the slow modes, 
and accordingly the slow variables of the two systems, are different, it is im- 
possible to say that the dynamic universality class of the two systems are 
identical. What this means is the breakdown of the prescription for choosing 
the slow variables just by gathering the order parameter and the conserved 
quantities, because the canonical momentum for the meson fleld is neither the 
order parameter nor the conserved quantity. In the chiral phase transition, the 
order parameter and the conserved quantity does not furnish us with the full 
slow variables required to describe the slow mode properly. We have danger- 
ous potentiality that the dynamic universality class given by Hohenberg and 
Halperin's classiflcation method is divided into flner classes, if we consider the 
slow modes and their appropriate slow variables in that dynamic universality 
class. 

The necessity of the canonical momentum for the slow variable is itself noth- 
ing strange. When we return to the microscopic eqution of motion, le., the 
Hciscnbcrg equation for the meson fleld, we can immediately realize that the 
canonical momentum in addition to the meson fleld is necessary to describe 
the meson dynamics. Moreover we have mentioned that the canonical mo- 
mentum as a slow variable is not restricted to the chiral phase transition. In 
the suprefluid and the ferromagnet, the canonical momentum plays the role 
of the slow variable in order to describe the second sound and the spin wave 
respectively. In the two systems, the canonical momentum is just an order 
parameter or a conserved quantity, and can become the slow variable just by 
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the old prescription. In the chiral phase transition, on the other, the canonical 
momentum is not the order parameter or the conserved quantity. This is just 
the only new feature of the chiral phase transition. 

The above observation has inspired the necessity of reanalysis of the chiral 
phase transition. We have employed the mode coupling theory for that pur- 
pose. 

Firstly, we have applied the theory to the 0(2) linear sigma model and found 
that the meson mode appears desirably from the meson field and the canonical 
momentum. We have also clarified the other slow modes than the meson mode: 
The chiral charge gives a diffusive mode in the disordered phase. In the ordered 
phase, it couples with the pion field, which results in lifting the degeneracy 
of the pion mode with the sigma mode. The energy and the longitudinal 
momentum are combined to give the energy wave mode, which is a propagating 
mode and a correspondent to a sound wave in a fluid. The transverse momenta 
give diffusive modes. 

We have then examined the 0(4) linear sigma model in the mode coupling 
theory in order to derive the kinetic equation and calculate the dynamic crit- 
ical exponents in the disordered phase. We have found the exponent for the 
meson mode to be z,^ = 1 — r]/2 = 0.98, which is to be compared with the 
value oi d/2 obtained assuming that the chiral phase transition belongs to the 
dynamic universality class of the antiferromagnet. The different dynamic crit- 
ical exponents show explicitly the different dynamic behavior. We have also 
calculated the dynamic critical exponents of the modes other than the meson 
mode. The results are given in Eqs. (6.54), (6.60) and (6.65). We have noted 
that the largest exponent is given by the transverse momentum fluctuation, 
which determines the relaxation time of the whole system in the chiral phase 
transition. 

We have succeeded in describing the meson mode by including the canonical 
momentum into the member of the slow variables. We were able to flnd, even 
accidentally, the needed slow variable, that is, the canonical momentum by 
returning to the microscopic Heisenberg equation. However, we must admit 
that the argument does not go beyond a heuristic one. We do not have any 
deflnite prescription to determine the slow variables. It will be our future 
problem to flnd out the way how to determine the slow variables uniquely, if 
it exists. Once it is found, it means that we have obtained a new classiflcation 
method for the dynamic universality class, which completes that of Hohenberg 
and Halperin. 

We note that a "particle mode" like the meson mode is inherent in the rela- 
tivistic system because it accompanies the anti-particle mode as a partner. As 
we have seen, for the description of the "particle mode," we need the canoni- 



46 



cal momentum that is neither the order parameter nor the conserved quantity. 
Thus the relativistic system that involves the "particle mode" as a slow mode 
should be a quite novel critical point that is not classified into any dynamic 
universality classes of the non-relativistic systems considered by Hohenberg 
and Halperin in Ref. [14]. One of such a system would be the critical end 
point (CEP) in the QCD phase diagram [29,30,31,32,33,34,35]. Although the 
dynamic universality class of the CEP is discussed in Ref. [35], the argument is 
not sufficient in respect that only the order parameter and the conserved quan- 
tities are compared and the slow modes are not taken into account explicitly. 
Other new and intriguing systems would be the tricritical point, the critical 
point associated with the color-superconductivity, the confinement transition, 
the electro-weak transition and so on. The analysis of those systems including 
CEP should also be what we should do in the future. 
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A Commutation relations of the slow variables 



In this appendix, we list up the commutation relations among the slow vari- 
ables. 
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2to 



aa' '->qq' 



+ 
+ 



to-, 



q-q' 



(A.l) 
(A.2) 

(A.3) ■ 



1 1 



t 

aq'—q 



A 



si- 



-g',fei+-+fc4'^-9,fci 



9 ki-h * V'/ci "fe4 

(A.4) 



'<l>LqA<f>l-q')i 
0i_„-, (<3a'/3'g')^ 



■aq-q' -r ^'^aq' -q 



+ 



r, ^aa' ^qq' 
ZUq' 

- [<i>a-q, 



(A.5) 
(A.6) 



Q-Q' I _ t 

UJr' 



(A.7) 
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1 



1 ^ / t 

1 



2 V V^^/2cJ^ y'cJ^u^ 

fei---A;4 



(A.8) 



4'a—qj Pqi 



daq-q' + O'aq'-q 



q~q' if f 

1 [^aq-q' " ^ag'-g' 



Qlq^ (Qlq')\ 



Qaq^ iQbq' 
[Qaq^ iQbq' 
[Qal3q, [Eq") 



i 

i 

"7v 



^abcQca—< 



cq-q' 



^abcQ, 



cq-q' 



^abcQcn- 



V 

cq-q' 
1 ^ 1 



(A.9) 
(A.IO) 

(A.ll) 

(A.12) 



= - le. 



q-q'-k 



X 



{q-k)-{q-^ -k) (a^g + a^_-) (a^-_-,r + a 



QaPq, (P^)^ 

E^, {E, 

Eq, (P^ 



(fQa/iq-q' 

{q + q'yPl_^ 



Sq-q'k ^ '^Sk+q'-q 

(A.13) 
(A.14) 

(A.15) 



1 11 



^ ^ak+q' 



+ {q-(^ -k) 



1 



A ^ 1 



fe: 



fel '^ifc4 



4"q-g',fci+-+fc4 



A;i---A;4 



«3 a-ksJ V a-fe4/ 



PhiPP'] =y{^PP^ + <^'P^- 



(A.16) 
(A.17) 
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B Dissipation terms in the relativistic hydrodynamics 



The dissipation terms for and Pi are read from the relativistic hydro- 
dynamic equation [36]. We will perform it in this appendix. The energy- 
momentum tensor for the imperfect fluid is given by 

T"'^ = To"'^ + AT"^, (B.l) 
where 

To"^ = prj^/^ + {p + p)U''U^, (B.2) 

AT""^ = -r]H'^^ H^^W^s - A (//°^C/^ + H^^V) - (H""^-^. (B.3) 

The Tq^ is for the perfect fluid and AT°^^ takes care of the dissipation effects. 

The p and p are the pressure and the proper energy density. The coefficients 
A, 7], ( are the heat conductivity, the shear viscosity and the bulk viscosity, 
respectively. The metric tensor r)°'^ are taken to be r)°'^ = diag(— 1, 1, 1, 1). 
The four velocity C/" = — v^, v/\/l — v"^) has the property U"Ua — — 1 

in this metric. The Wap, Qa, H^/s are defined by 

HaP = Tlap + UaUp. (B.6) 

We define the energy density and the momentum density as 

E(f,t)=Tr = ^±^, (B.7) 

i. — 

We can derive the relation between Lf^ and P*. ^Prom Eq. (B.8), we see 

P'^{p + p)U^U\ (B.9) 
Moreover using C/"C/q = —1, we find 



in the approximation of the small momentum fluctuation P* <^ 1. 
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Now consider the hydrodynamic equation. It is given by the energy-momentum 
conservation, 



= 



from which we find 



dt 



dx^ 



dx^ 



dx^ ' 
dAT'^ 

dxP ' 



(B.ll) 

(B.12) 
(B.13) 



The second terms in the right hand sides are the dissipation terms, which we 
will compute in the following. The AT°'^ are calculated to be 



BT BT 



dx^ 



dx 



dXr 



dxp dxa 



dx-y 



dx^ 
dx^ 



dx^ ' 



(B.14) 



Thus the dissipation term for the energy density is 



dx^ 



dx^ 



+ 2W 



dxo 



+ TW 



+C (-i + c/°f/^ 



dx"* J 



dx-y 



d 



dx^ 



+r] 



A \2U^WU"' 



+ 



dT 

dx"/ 



dT 

dxj 



dx, 







dxi 



dxj dxo 



3 dx^ dx-^ 



(B.15) 



Here the temperature T is related to E through the thermodynamic relation 



5T = -5E 



(B.16) 



with C being the specific heat. The and are translated to by Eqs. 
(B.9) and (B.IO). Moreover following the non-relativistic case, we take only 
the terms up to the third order of the assumed small quantities, that is, the 
frequency and wavenumber as well as the fiuctuations E and P\ Since each 
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term in Eq. (B.15) involves two time and space derivatives or the frequency 
and wavenumber, we have only to take the first order of E and P\ Thus we 
find 



= A 
^ A 



d 

d 

dx^ \ cdx, 



1 dE ^^d 



C dxj 
' 1 dE 



.pj 



dt \p + p 
1 (9P^~ 



Po + Po dt 



(B.17) 



where in the last step, we replaced T, p, p with their equilibrium values Tq, 
Po, Po- 



Similarly, the dissipation term for the momentum density becomes 



1 dE 



+ To- 



Qpi 



+ v 

+ C 



Cdxi Po + Po dt 
1 d 



Po + Po dx^ \ dxj 

1 d dP^ 
Po + Po dxi dx^ ' 



dP' dP^ 

+ 



dxi 



3 ' dx^ 



(B.18) 



The hydrodynamic equation with the dissipation term is thus 
dE___dP^ _9_ A as To dP^ \ 

'dt~~'d^^ d^ [cd^j^ Po + po dt ) ' ^ ' 

Qpi Qjdj Q ( I QE To dP'\ 
+ A— * 



dt dx^ dt \C dxi po+Po dt J 

1 d fdP' dP^ 2 ,-,9P'\ ,1 d dP 
+ V — : — TT- + ^ - 7;V''^ + C 



Po + Po dx^ \ dxj dxi 3 dx'' J Po+ Po dxi dx^ 

(B.20) 

We note that there appear the relativistic corrections with the time derivative. 

If those terms are dropped, the dissipation terms reduce to the non-relativistic 
ones. Moreover we note that the thermal conductivity A enters into the equa- 
tion for as a relativistic effect. 

When moving to the Fourier space, we finally obtain 

dE, . . / 1 2p ^ To . dPl 



= - i?.-^^ + A -^q'E^-r -r^iq,^ , (B.21) 



dt • \ " • po + Po ^' dt 

..rr^ii^J^.^^dE. To ^P^ 



dt J y (J jj^ jj^2 

+ % in'P'^ + k'P^ - y^mP^) + ^^iq'iqiPi^, (B.22) 

Po+Po V ^ '3 V po+Po ^ 
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from which we find the dissipation terms in Eqs. (5.48), (5.55). 



C Scaling properties of the commutators of the slow variables 



In this appendix, we examine the scahng properties of the commutators of 
the fluctuating amphtudes, for {(!),(})'', Q,E,P''} in the chiral phase 

transition. The results provide us with the exponents Xji, which are used to 
find the scaling properties of the commutators of the reduced amplitudes. 

The commutation relation is given in Eq.(A.l). 



(t>aq, (0qy)^ 



1 



(C.l) 



where we have used cUq = {2x<i}q) ^^'^ ~ C 
This is zero. 

• [(Paq, {Qa'fS'q'Y] 

The commutation relation is given in Eq. (A. 3). Noting that 



(C.2) 



we find 



<Paq, {Qa'P'q'y 



—1 



2VV 



(C.3) 



[(t>aqAEq'y] 

Using the approximation [Ag, Sj^,] ~ [A^_^-^, B^, _^ with |Z| <S 1, 



we find 



^aq—q' -i 



H 



namely 



(C.4) 



(C.5) 
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The scaling property of can be found by considering 

Denoting the exponent of as x^, we find 
which gives x^ = d/2 + (3/v. Thus scales as 
leading to 



(C.6) 

(C.7) 

(C.8) 
(C.9) 



From Eq. (A.5), 



This gives the scaling properties, 



(C.IO) 
(C.ll) 



for the longitudinal and transverse components of the momentum. 

[QalBq-, {Qa'jB'q'V] 

The commutation relation Eq. (A. 10) immediately gives 



Q^q^ {Qbq' 



^ahcQcq-q' ~ ? 2 ■ 



(C.12) 



As before, the commutator is approximated as 



[Qal3q,E_^] ~ [Qapq-q'.E^] 



Qafiq—q'-i ^j^"^ 



Qa, 



that is. 



[Qa/J?, E. 



(C.13) 
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Consider 



■k-k' 



^'^'^ (~7yQk+k'Q-k-k' + TTy^k'Q-k' 



(C.14) 

This gives 

lid.-XQ-2xQ ^ ^-f j^d-2xQ (C.15) 

to find Xq = d where Xq denotes the exponent for Q. The scahng property 
of Q is thus 

from which we obtain 

[Q„^,-, E_,.] = -^Q^p {q- ^, V) ~ r'^. 

The commutation relation Eq. (A. 14) gives 



The commutation relation Eq. (A. 15) gives 



The same procedure as before leads us to 



pi 

V '-'^ 



that is, 



—1 



-k-k'/ ^ V k+k'^ -k-k' 



(C.16) 



(C.17) 



(C.18) 



(C.19) 



(C.20) 



Consider 

H4P-k-k)-k^(\pi^4P-k-k] 
-^{k^{pu,pL 

+ >^{n'P-k) + (k + J^y{PiPLk)}, (C.21) 
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which gives the exponent for P = d + 1. Thus 



and we find 



9' ^ 9' 

The commutation relation Eq. (A. 17) gives 
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